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ABSTRACT 



Aims. We have studied the dependence of the profiles of molecular abundances and line emission on the accretion flow 
in the hot lOOK) inner region of protoplanetary disks. 

Methods. The gas-phase reactions initiated by evaporation of the ice mantle on dust grains are calculated along the 
accretion flow. We focus on methanol, a molecule that is formed predominantly through the evaporation of warm ice 
mantles, to show how the abundance profile and line emission depend on the accretion flow. 

Results. Our results show that some evaporated molecules keep high abundances only when the accretion velocity is 
large enough, and that methanol could be useful as a diagnostic of the accretion flow by means of ALMA observations 
at the disk radius of ^ lOAU. 
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1. Introduction 

Recent observations have detected a variety of molecu- 
lar lines towards disks around T Tauri stars (e.g., Dutrey 
et al. 1997; Thi et al. 2004; Qi et al. 2008). Existing 
millimetre/sub-millimetre observations, which have rela- 
tively low spatial resolution, can trace only the outer re- 
gion 50AU) of the disk, while near-infrared observa- 
tions with high sensitivity or high spectral resolution probe 
molecular lines from the planet-forming region in the disks 
(Lahuis et al. 2006; Gibb et al. 2007; Carr & Najita 2008). 
The forthcoming Atacama Large Millimeter/sub-millimeter 
Array (ALMA), with both high sensitivity and high spatial 
resolution, will make it possible to observe various molecu- 
lar lines from the inner regions of the disks. 

Although many models on the chemical structure of 
young circumstellar disks focus on the outer disk (e.g., 
Aikawa et al. 2002; Willacy 2007), the chemistry in the 
inner disks has also been studied (e.g., Markwick et al. 
2002). In cold pre-stellar cores, observations suggest that 
many gas-phase molecules, including CO, are frozen onto 
dust grains (e.g., Caselli et al. 1999). Young stars and cir- 
cumstellar disks are thought to be formed as a result of the 
collapse of such molecular cloud cores. In addition, observa- 
tions of molecular lines towards disks suggest that the gas- 
phase molecules are frozen onto grains near the midplane 
of cold outer disks (e.g., Dutrey et al. 1997). However, in 
the inner region of the disk the dust temperature is high 
due to the irradiation from the central star, and the ice 
mantles are expected to evaporate into the gas by analogy 
with the molecules observed towards star-forming cores, so- 
called hot molecular cores and hot corinos (e.g., Millar 1993; 
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Ceccarelli et al. 2007), so that the chemistry in the inner 
disk will be characterized by this ice mantle evaporation. 

Accretion flow towards a central star is important for 
planet formation in the disk because it is one of the possible 
mechanisms for the dispersal of the gas, (e.g., HoUenbach 
et al. 2000), while the amount of the gas in the disk influ- 
ences the coagulation and settling of dust grains, which lead 
to planetesimal formation, and also controls gaseous planet 
formation as well as the orbital motion of the planets (e.g., 
Nakagawa et al. 1986; Kominami & Ida 2002; Papaloizou et 
al. 2007). Meanwhile theory has suggested that the planet 
forming region in the disks could be magnetorotationally 
stable, that is, there could be no driving source of the ac- 
cretion flow (e.g., Sano et al. 2000). Therefore, it will be 
useful if we can identify observational evidence for the ac- 
cretion flow in the inner disk in order to understand the 
planet formation processes. 

If the gas and dust in the disks are accreting towards 
the central star, the ice mantle molecules, which are frozen 
onto dust grains in the cold outer disk, will evaporate when 
they are transported into the hot inner disk. The evap- 
orated molecules will be destroyed by chemical reactions 
while they are migrating inwards along the accretion flow. 
Thus, the distributions of the molecules will trace the ac- 
cretion velocity, depending on whether the timescale of the 
chemical reactions is shorter or longer than the accretion 
time. Hence, it will be possible to check the existence of 
accretion flow from observations of the molecular distribu- 
tions in the inner disk. 

In this work we study the effect of the accretion flow on 
the chemical structure of the inner region of young disks by 
modelling the chemistry, initiated by the ice evaporation, 
along the accretion flow. And we examine the effect on the 
observational properties, especially focusing on a transition 
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Fig. 1. Contour plots of the gas temperature {solid lines) 
and density (dotted lines) distributions in the z/x vs. x 
plane for the fiducial model with M = 1.0 x 10^^ Af© yr~^. 

line of methanol, which is an abundant ice molecule in many 
high and low mass star forming cores (e.g., Gibb et al. 2004; 
Boogert et al. 2008). Methanol, which is difficult to form 
efficiently in the gas phase, is also known to be more abun- 
dant in the gas-phase around high and low mass young stars 
as well as comets (e.g., Macdonald et al. 1996; Schoier et al. 
2002; Bockelee-Morvan et al. 1991) by more than orders of 
magnitude compared with cold dark clouds (e.g., Ohishi et 
al. 1992). In addition, laboratory experiments have shown 
that methanol can be formed in ice mantle with reason- 
ably short timescale as a result of hydrogenation of carbon 
monoxide (e.g., Watanabe et al. 2006). Therefore, methanol 
is thought one of the most probable components of grain 
mantle molecules. 

In the following section we introduce the physical and 
chemical models of the disk. We present the resulting molec- 
ular abundance profiles of the inner disk in Sect. 3 and the 
line emission of methanol in Sect. 4, using the models with 
different accretion velocities. Influence of some assumptions 
on the results are discussed in Sect. 5, and the results are 
summarized in Sect. 6. 

2. Models 

We model an axisymmetric disk surrounding a central star 
with a mass of M* = 1.5Mq, a radius of i?, = 2i?0, and 
a temperature of T, = 6000K, whose luminosity (~ 5Lq) 
is relatively higher than the averaged luminosity of typical 
T Tauri stars (- ILq) (e.g., D'Alessio et al. 2005). It is 
observationally known that many T Tauri stars emit strong 
X-ray and UV radiation. We adopt a model with black body 
(T, = 6000K) plus thermal bremsstrahlung (Tbr = 2.5 x 
lO^K) emission for UV radiation (Lpuv — 7 x lO'^^ergs s~^; 
see Nomura & Millar 2005, hereafter NM05) and a simple 
thermal bremsstrahlung emission model with Tx = IkeV 
and Lx = 10^" ergs s~^ for X-rays (Glassgold et al. 2004; 
Garmire et al. 2000). 

The temperature and density distributions of gas and 
dust in the disk are obtained self-consistently as described 
in previous papers; see NM05 and Nomura et al. (2007) for 
details. Dust temperature is obtained by the iterative radia- 
tive transfer calculation (see Nomura 2002), where an ini- 
tial dust temperature profile is calculated using the variable 



Table 1. Initial fractional abundances with respect to total 
hydrogen nuclei. 



Species Abundance 


Species 


Abundance 


H+ 


1.0 X 10' 


ii 


CO 


1.3 X 10" 


4 


He+ 


2.5 X 10" 


12 


CO2 


3.0 X 10" 


6 




1.0 X 10' 


-9 


H2CO 


2.0 X 10" 


6 


Fe+ 


2.4 X 10- 


-8 


CH3OH 


2.0 X 10" 


7 


He 


1.0 X 10' 


-1 


C2H5OH5.O X 10" 


9 


S 


5.0 X 10" 


-9 


O2 


1.0 X 10" 


6 


Si 


3.6 X 10- 


-8 


H2O 


2.8 X 10" 


4 


C2H2 


5.0 X 10^ 


-7 


N2 


3.7 X 10" 


5 


CH4 


2.0 X 10" 


-7 


NH3 


6.0 X 10" 


7 


C2H4 


5.0 X 10" 


-9 


H2S 


1.0 X 10" 


7 


C2H6 


5.0 X 10" 


-9 


ocs 


5.0 X 10" 


8 



Eddington factor method (DuUemond et al. 2002) in this 
work. The disk is assumed to accrete towards the central 
star with a constant mass accretion rate, and the a-viscous 
model with a = 0.01 is adopted. The dust model which re- 
produces the observational extinction curve of dense clouds 
(see NM05) is adopted here. The resulting gas temperature 
and density profiles for a fiducial model with the mass ac- 
cretion rate of i\f = 1.0 x 10^^ Mq yr^^ are plotted in 
Figure [1] 

For the chemical model we calculate the time-dependent 
gas-phase reactions along stream lines of steady accretion 
flow as 

'gf^ = '^'i^^ + ktjkn^nk, (1) 

where Ui is the number density of the species, i, and the 
right hand side of the equation describes the formation 
and destruction of the species, i, due to cosmic-ray /photo- 
chemistry (the first term) and two-body reactions (the sec- 
ond term). Three-body reactions are not included, since 
they will not affect the result very much in the gas temper- 
ature range we treat in this work (e.g., Willacy et al. 1998). 
We assume that the stream lines, I, are simply described as 
z = sH, where H — Csq/^k is the disk scale height, and 
Cso the sound speed at the disk midplane, fix the Keplerian 
frequency. We set 60 grids for < s < 3. The equations 
are solved together with the continuity equation along the 
stream lines. The accretion velocity, Wacc, is simply given by 
Vacc = M/{2tt'Ex), where M is the mass accretion rate, E 
the surface density of the disk, and x the distance from the 
central star. Turbulent mixing is not taken into account in 
this paper (see discussions in Sect. 5). 

The chemical network consists of 208 species and is 
connected by 2830 reactions, in which the reaction rate 
coefficients are taken from the UMIST RATE06 database 
(http://wwwjadfa.net/) (Woodall et al. 2007). The X-ray 
ionisation is simply modelled by analogy with the cosmic- 
ray ionisation (see e.g., Aikawa et al. 1999 for a more de- 
tailed model). X-rays will ionise the gas and induce photore- 
actions, similar to cosmic-rays; we adopt the rates enhanced 
by Cjf/CcR for the reactions, where the X-ray ionisation 
rate, (x, is calculated using the X-ray flux at each position 
in the disk (see Nomura et al. 2007). The H2 cosmic-ray 
ionisation rate is set to be (qr = 1-3 x 10"^'' s~^, and 
the attenuation is neglected since it is not significant ex- 
cept very close to the disk midplane at the disk radii of 
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X ~ lAU in this model. As the initial condition of the 
calculations, the fractional abundances of the species are 
taken from Nomura & Millar (2004; hereafter NM04), in 
which the initial ice mantle composition is fixed so as to be 
consistent with the infrared absorption features of ices ob- 
served towards young stellar objects as well as the line emis- 
sion of gas-phase molecules observed towards a hot core, 
G34.3-h0.15 (see Table 1). 

Since we particularly focus on the observational prop- 
erties of a methanol line in Sect. 4, we simply start the 
calculations where most methanol evaporates into the gas, 
that is, where the timescale of the thermal evaporation 
of methanol, Tovap, becomes shorter than its adsorption 
timescale onto dust grains. Tads (see Fig. [5]). The evapo- 
ration time is given by Tovap — eyij>[Ei, / kTd) , where 
vq is the vibrational frequency of the adsorbed species on 
the grain surface. The timescale is very sensitive to the 
binding energy of the molecules onto the dust grains Ef,, 
and the thermal energy of the grains, kTd- The binding en- 
ergy of methanol of Eh/k — 4235K is used in this work, 
which is experimentally measured for pure methanol ice 
(Sandford & Allamandola 1993). We note that the bind- 
ing energy of methanol on methanol is equal to that of 
methanol on water ice within the error bars (W.A.Brown, 
private communication). The adsorption time is given by 
Tads — {Sira^dgnvth)^^ , where S is the sticking probabil- 
ity, TTo^dgn the grain surface area per unit volume (dg is 
the abundance of the dust grains relative to total hydrogen 
nuclei), and Vth the thermal velocity of the molecule (e.g., 
Tielens & Allamandola 1987). Here we use S = 0.3 and 
the averaged grain surface area per total hydrogen nuclei 
of < dgO^ >= 2.2 X 10-22 cm which is similar to that of 
the MRN dust model by Mathis et al. (1977) (see NM04 
for details of the parameters; see also discussions in Sect. 5 
for effect of dust evolution). In this model all the ice man- 
tle molecules evaporate at the same dust temperature of 

~ 85 — lOOK, depending on the local density. The influ- 
ence of this simple treatment of the ice mantle evaporation 
is briefly discussed in Sect. 5. 
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Fig. 2. Evolution of molecular abundances relative to total 
hydrogen nuclei along the accretion flow, z 1.2H, for a 
model with the mass accretion rate of M = 1.0 x lO~^Af0 
yr"^. The figure at the bottom shows the gas density (n) 
and temperature (T) profiles and the distance from the 
central star (x) . The molecules evaporated from ice mantles 
arc destroyed by chemical reactions around the timescale of 
- lOV- 



3. Resulting molecular distributions 

We calculate molecular abundances using the obtained den- 
sity and temperature profiles. Figure [2] shows the resulting 
evolution of abundances of some molecules relative to total 
hydrogen nuclei along the stream line of z = 1.2H. The 
results for a fiducial model with the mass accretion rate 
of M = 1.0 X 1O-^M0 yr"^ are plotted, in which the gas 
density and temperature range over n ^ 3 x 10^^ — 1 x 10^'^ 
cm"^ and T ~ 100 — 240K, respectively, along the flow. 
The gas density (n) and temperature (T) profiles as well 
as the distance from the central star (x) are also plotted in 
the figure at the bottom. 

Figure [2] shows that the parent molecules which evap- 
orate from grain surface, such as NH3, H2S, C2H2 and 
CH3OH in this model, are relatively stable, but are eventu- 
ally destroyed by reactions with ionised species or atomic 
hydrogen to produce daughter species, such as HCN, SO2, 
and CH3OCH3, around a timescale of Tchcm ^ lO^yr. The 
parent and a part of daughter species end up as CO and 
small hydro-carbon molecules. The timescales for the de- 
struction of these parent molecules are nearly independent 
of the total density, n. This is because the number densities 



of ionised species and atomic hydrogen, which are formed 
via cosmic-ray ionisation and induced photoreactions, re- 
spectively, are not very sensitive to n. OCS is destroyed 
on a longer timescale than other parent molecules since its 
main reactant is H;^, while many of the other molecules 
are destroyed mainly by H3O"'" which is more abundant 
than H;^. OCS does not react with H3O+ because its pro- 
ton affinity is smaller than that of water. H2O and CH4 are 
hardly altered since they are produced efficiently in the gas- 
phase, especially at high temperature in the case of H2O, 

which is formed via endothermic reactions, O OH 
H2O. Methane is formed via the destruction of methanol 
and the formation processes are less temperature depen- 
dent. The abundances of NH3 and HCN increase around 
t ~ lO^yr in Figure [2] owing to endothermic reactions and 
high temperature in the inner disk. The detailed gas-phase 
chemistry induced by the evaporated molecules is similar to 
that occurring in hot molecular cores, which is discussed in 
NM04. Along z = 1.2H molecules are not affected by UV 
and X-ray irradiation from the central star in this model 
(see also Sect. |4]). We note that the timescales for destruc- 
tion of these parent molecules by chemical reactions have 
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uncertainties of about an order of magnitude due to the 
uncertainties in the rate coefficients of the reactions (e.g., 
Wakalem et al. 2005). 

In Figure [31 we plot the abundance profiles of some 
species relative to total hydrogen nuclei as a function of 
the disk radius along the stream line of 2; = 1.2H for mod- 
els with mass accretion rates of M = (a) 1.0 x 10~^ and (b) 
5.0 X 10~^Mq yr~^. In the steady accretion flow, a location 
of a fluid particle in the disk, x, is simply related to a time 
after the ice evaporation, t, as t = {vaccix')}~^dx' 

(where a^evap is the evaporation radius), as shown in the 
bottom of Figured! The radial dependence of the accretion 
velocity is weak, and Wacc ~ 40 cm s~^ for the fiducial model 
with M = 1.0 X 10~®Mq yr~^. In this work the accretion 
velocity is set to be proportional to the mass accretion rate 
('^acc oc M) by artificially fixing the surface density of the 
disk to that of the fiducial model for the purpose of clear 
understanding of the effect of accretion flow on the chemical 
structure in the disk. 

Figure [3] shows that the molecular abundance profiles 
in the inner disk change dramatically depending on the ve- 
locity of the accretion flow. When the accretion velocity 
is small and the accretion time, Tacc ^ x/v^cc, is longer 
than the timescale of the chemical reactions which de- 
stroy the evaporated molecules, Tchcm, the parent molecules 
are transported only a small distance along the accretion 
flow before they are destroyed by the chemical reactions. 
Therefore, some parent and daughter species are abun- 
dant only in the region where parent molecules are evap- 
orated into the gas. On the other hand, if the accretion 
velocity is high and the accretion time is shorter than the 
chemical timescale, the parent molecules are transported 
to the central star before they are destroyed, and high 
abundances occur even close to the star. Some daughter 
species, such as HCN, become abundant gradually as the 
gas is transported inwards in this case. We note that the 
abundances of some daughter species are very sensitive to 
the initial fractional abundances; for example, if atomic 
nitrogen is abundant initially, the maximum value of the 
HCN abundance increases by about two orders of magni- 
tude compared with the model we used here. Along the 
stream line of z = 1.2H, C2II2 is destroyed by an endother- 
mic reaction, C2H2+O— s-C0-|-CIl2. Thus, the C2H2 abun- 
dance decreases, independent of the accretion rate, around 
a; ~ 4 — 5AU, where the gas temperature becomes high 
enough for the above reaction to proceed. Along stream 
lines in upper layers of the disk slightly above z = 1.2H, 
where the gas density is relatively low, C2H2 is destroyed by 
reactions with molecular ions around a timescale of lO^yr, 
similar to other parent molecules. 

4. Methanol line emission 

Using the gas density and temperature as well as the 
molecular abundance profiles obtained in Sect. 3, we cal- 
culate the brightness temperature of 145GIIz Jk = 3o — 2o 
transition of methanol. Figure |4] shows the resulting 
brightness temperature of the line along a single line of 
sight as a function of the disk radius for models with 
the mass accretion rates of M = 1.0 x 10~^ {solid line), 
1.0 X 10~^ (dotted line), and 1.0 x 10~^ {dot-dashed line) 
Mq yr~^. We assume that the molecules are in the local 
thermal equilibrium. The assumption is reasonable since 
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Fig. 4. Brightness temperature of the 145GIIz Jj^ — 3o — 2o 
line of methanol calculated along a single line of sight as a 
function of the disk radius for models with mass accretion 
rates of M 1.0 x lO"'^ {solid line), 1.0 x 10"^ {dotted 
line), and 1.0 x 10"^ {dot-dashed line) Mq yr~^ The ALMA 
detection limit for the beam size of ~ 14AU is also plotted 
with a thin solid line, showing that the line emission will 
be detectable in the inner disk when the accretion rate and 
the accretion velocity are high. 



the line is emitted from regions where the density is 
much higher than the critical density for this transition, 
Ucr 2 x lO^cm"^ (e.g., Pavlychenkov et al. 2007). 
The disk is assumed to be face-on to an observer. The 
molecular data is taken from Leiden Atomic an d Molecular 
Database ( ht t p : / / www . st r w . leidenuniv . nl/ ~ | moldata/) 
(Schoier et al. 2005), ^WKflT part of the RATRAN code 
(http://www.sron.rug.nl/"' vdtak/ratran/ratran.html) 
(Hogerheijde & van der Tak 2000) is used for the line 
radiative transfer calculations (see NM05 for the detailed 
equations used in the calculations). 

The figure shows that the brightness temperature of the 
line, Tb, is high all over the inner disk when the accre- 
tion velocity is high, while Tb becomes high just around 
the region where the methanol is evaporated into the gas 
when the accretion velocity is low, since Tb traces the abun- 
dance of methanol (see Sect. 3). Now, methanol is mainly 
destroyed by chemical reactions with II3O+, on a timescale 
around Tchcm ~ (fc"^H30+)^^ lO^yr. Here the rate coef- 
ficient is fc ~ 3 X 10~^cm^s~^ and the number density of 
H3O+ is nH30+ ~ 10"'^ cm"'^, which is almost independent 
of the total number density, n, since the formation process 
of H3O+ is related to cosmic-ray ionisation. The timescale 
of methanol destruction is slightly affected by the abun- 
dance of NH3 which helps recycle methanol by reacting 
with CH3OH+ (e.g., Rodgers & Charnley 2001). When the 
accretion rate is high, Tb is high even in slightly outer re- 
gions {x ~7-9 AU) because the dust temperature is higher 
due to the viscous heating. 

At the disk surface methanol is photodissociated due 
to UV and X-ray irradiation from the central star. The 
timescales of the dissociation are roughly given by ruv ~ 
fcyv ~ SOGpuv yr and tx ~ ~ lO'^Cx^ respec- 
tively, where Gpuv is the FUV radiation field normalized by 
1.6x 10"'^ ergs cm~^ and is the total hydrogen ionisa- 
tion rate by X-rays. Therefore, the photodissociation dom- 
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Fig. 3. Abundance profiles of some species relative to total hydrogen nuclei as a function of the disk radius along the 
stream line of z = 1.2H for models with the mass accretion rates of M = (a) 1.0 x 10^^ and (b) 5.0 x 10~^Mq yr~^. 
Some molecules are abundant just around the region of the ice mantle evaporation when M is low and the accretion time 
(■Tacc) is longer than the timcscalc of the chemical reactions which destroy the parent molecules, while they are abundant 
even close to the central star when M is high and race is short enough. 



inates the destruction of methanol where Gpuv ^3x10 ^ 
or Cx lO"^'"^ s^^. In Figure [5] the contour lines of 
Gpuv = 10"^ and lO^^" {thin solid lines), Cx = 10"^^ 
and 10^^^ s^^ (dashed lines) are plotted in z vs x plane of 
the disk together with the line of Tovap = Tads for methanol 
{thick solid line), which corresponds to ~ 85 — lOOK, 
and z = H {dot-dashed line) for the fiducial model. Figure[5] 
shows that the line brightness temperature at a; > 6— 7.5AU 
in Figure[3]is low because the methanol can evaporate from 
dust grains only in the surface layer of the disk where it 
is dissociated quickly by strong X-ray and UV irradiation 
from the central star. 

In Figure [4] we also plot with a thin solid line 
the ALMA detection limit for the methanol line, cal- 
culated using the ESO AL MA Sensitivity Calculat or 
(http: / / www.eso.org/sci/facilities / alma/ observing/tools 
/etc/index. html). The detection limit of ^ 17K is achieved 
when we observe for 600 sec by using 50 antennas of 
the 12m Array, assuming a velocity resolution of 1 km 
s~^, and a beam size of 0.1 arcsec, which corresponds to 
~14AU when the disk is located at the distance of nearby 
star-forming regions, ^ 140 pc. The 145GIIz Jk = 3o — 2o 
transition of methanol has not been detected towards 
DM Tau (Dutrey 2001), and other methanol line searches 
towards some T Tauri and Herbig Ae stars (Thi et al. 
2004; Semenov et al. 2005), have also been unsuccessful. 
Our results suggest that when the accretion velocity at 
X < lOAU is high enough, ALMA can detect the methanol 
line, taking advantage of high spatial resolution to avoid 
beam dilution. When the accretion velocity is low enough, 
line detection will be difficult even with ALMA. Thus, the 
observations of methanol by ALMA could be a useful tool 
to diagnose the existence of accretion flows in the inner 
disks. 

5. Discussions 

In this section, we discuss some possible effects of the as- 
sumptions on our results. 

First, we discuss the influences of the simple treatment 
of the evaporation process we have taken in this work. 
Although we simply assume that all the species evaporate 
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Fig. 5. Line of 

'^cvap = Tads inside which the icy methanol 
evaporates into the gas {thick solid line), line oi z = H {dot- 
dashed line), and contour lines of Gpuv = 10^'^ and 10^^" 
{thin solid lines), Cx = 10^^'^ and 10^^^ s^^ {dashed lines) 
for the fiducial model with M = 1.0 x lO^^M© yr~^ The 
evaporated methanol is quickly dissociated due to strong 
UV photons or X-rays at a; > TAU. 



at the same dust temperature, the evaporation of ice man- 
tle molecules may depend on the binding energy of each 
molecule and the morphology of water ice (e.g., CoUings et 
al. 2004). Therefore, different parent molecules may evapo- 
rate at different positions (see e.g., Markwick et al. 2002), 
which will affect the spatial distributions of molecular abun- 
dances. In particular, CO and N2 are likely to evaporate at 
temperatures much lower than 85K, particularly if they are 
abundant in the surface layers of the ice. However, given the 
chemical stability of these species, they are unlikely to affect 
the overall abundance of a molecule such as methanol. So, 
the timescales for the destruction of the parent molecules 
by chemical reactions, Tchcm, analysed in Sect. 3 will not be 
very sensitive to the treatment of the ice mantle evapora- 
tion. In addition to the evaporation process, surface chem- 
istry involving heavy radicals on warm dust grains may 
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affect molecular abundances in both solid- and gas-phase 
(e.g., Garrod & Herbst 2006; Aikawa et al. 2008). More 
detailed model including the detailed evaporation process 
as well as surface reactions on dust grains should be con- 
structed in future work. 

Next, if we take into account the effect of turbulent 
mixing in the vertical direction, the methanol abundance 
may increase near the line, Tevap = ^ads, (see Fig. [5]) when 
the diffusion time is shorter than the timescale of the pho- 
todissociation (e.g., Willacy et al. 2006). In the inner disk 
[x ^ 6AU) the methanol abundance may slightly decrease 
via the vertical mixing because of the photodissociation at 
the disk surface. The turbulent mixing in the radial direc- 
tion may enhance the gas-phase methanol abundance out- 
side the evaporation radius if the timescale of the turbulent 
diffusion is shorter than the timescale of the chemical reac- 
tions and that of the adsorption on dust grains. A detailed 
mixing model which treats angular momentum transfer will 
be needed in order to obtain the diffusion timescale, and 
advance further discussion. 

Finally, we note that the brightness temperature of the 
methanol line will reflect some conditions other than the ac- 
cretion velocity: the high-Te region becomes smaller if the 
central star is less luminous and the dust temperature in 
the disk is low. Also, Tb becomes low if the methanol abun- 
dance in the ice mantle is low. If the amount of small dust 
grains decreases due to the dust coagulation and settling in 
the disk, methanol is photodissociated by ineffectively at- 
tenuated FUV photons from the central star (e.g., Aikawa 
& Nomura 2006), and therefore Tb drops, although the 
evaporation radius moves slightly outwards owing to the 
decrease of the total grain surface area per unit volume, 
that is, the decrease of the adsorption rate (e.g., Aikawa 
2007). Thus, some causes other than low accretion velocity 
can lead to low methanol abundance and the non-detection 
of the line. However, if the observed methanol abundance 
is high in the inner disk, then the accretion velocity must 
be high in the planet-forming region. 

6. Conclusions 

We have calculated the gas-phase chemical reactions, initi- 
ated by the ice mantle evaporation, along the accretion flow 
in protoplanetary disks to show that the accretion velocity 
affects the abundance profiles of the evaporated molecules 
and their daughter species in the inner disk. We have ob- 
tained the brightness temperature profile of an emission 
line of methanol, and suggested that this could be useful 
as an indicator of accretion velocity in the planet-forming 
regions in disks using ALMA observations. 
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